Heat Capacity and Entanglement Measure in a simple two-qubit model 
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A simple two-qubit model showing Quantum Phase Transitions as a consequence of ground state 
level crossings is studied in detail. Using the Concurrence of the system as an entanglement measure 
and heat capacity as a marker of thermodynamical properties, an analytical expression giving the 
latter in terms of the former is obtained. A protocol allowing an experimental measure of entangle- 
ment is then presented and compared with a related proposal recently reported by Wiesniak, Vedral 
and Brukner 



PACS numbers: 

Introduction 

Entanglement as a strikingly non-classical feature 
has recently been linked to Quantum Phase Transitions 
(QPTs) in multipartite systems [T]-^. By definition, 
a quantum phase transition |3j is a sudden change 
of some physical properties a system exhibits at zero 
temperature, due to an abrupt modification of its ground 
state when a non-thermal parameter characterizing its 
Hamiltonian crosses a critical value [S] (known as quan- 
tum critical point or QCP). There exist two main classes 
of QPTs: the first type, known as First Order QPT 
(I-QPT), is a consequence of the occurrence of a level 
crossing [S] in the ground state energy in correspondence 
of the quantum critical point. This leads to a jump 
in the derivative of the energy of the system and to 
a discontinuity in the order parameter characterizing 
the transition. The second kind of QPTs, known as 
Second Order QPT (H-QPT) [7], can occur only in the 
thermodynamic limit when the dimension of the involved 
Hilbert space goes to infinity. In this case the quantum 
critical point is the limit of an avoided level crossing 
as the number of particles diverges. In both cases 
however the QCP manifests itself as a non-analytical 
point in the ground state energy. Consequences of such 
a non-analyticity, in some particular cases, might result 
in some peculiar properties of the system even at finite 
temperatures. 

A great effort has thus been devoted to the study 
of a particular kind of entanglement, namely the so- 
called thermal entanglement 8]. Its behavior, when 
non-thermal parameters are varied, is a signature of 
the physics underlying the transition itself [5]. It has 
been shown [10, that, under appropriate conditions, 
the occurrence of a quantum critical point against a 
non-thermal parameter k, at /t = k, is a necessary and 
sufficient condition for the existence of a discontinuity at 
K in some entanglement measures or in their derivatives. 
Moreover the renewed interest towards quantum thermo- 
dynamics in the last years [11] [T^ has spurred to search 
a possible link between the entanglement exhibited by 
a physical system in a thermal state and its thermal 
properties [I3HI5] . 

In view of these considerations it seems interesting to 



investigate the possibility of using a thermal quantity, 
like for example heat capacity, as a marker of a wide 
class of quantum phase transitions and to give a direct 
link between such a quantity and thermal entanglement 
in the state of the system under scrutiny. 
Some results in this direction have already been obtained 
[l6] but a full clear picture is not yet available. The 
aim of this work is to show that such an explicit link 
can be found, at least in a simple system of two spin-^ 
particles. 

This paper is organized as follows. In Section 1 we 
study a simple model of two qubits interacting with 
Heisenberg-like exchange potential and we define thermal 
quantities and entanglement properties for a thermal 
state of such a system. The main result of this paper 
is shown in Section 2, namely a direct link between 
a measurable thermal quantity and an entanglement 
measure, showing how this link can be employed to 
define an entanglement measuring protocol. Finally a 
comparison between our results and a similar (but not 
equivalent) result recently obtained by Wiesniak, Vedral 
and Brukner [15j is performed. 



I. HAMILTONIAN MODEL 

One of the most easily solvable spin models is the one 
described by a Heisenberg-like spin-spin interaction |17j . 
Its Hamiltonian characterizes many real physical situa- 
tions and it is perhaps the simplest Hamiltonian operator 
showing quantum critical points [TSl [TS] . 
Nevertheless, when the number of atoms grows an ex- 
act diagonalization of the corresponding Hamiltonian be- 
comes more and more difhcult, if possible at all. Even in 
a few cases where exact results are available in the ther- 
modynamic limit |20j . the dimension of involved Hilbert 
spaces is so high that analytical results are usually hard 
to be exploited in order to perform a detailed study j21j . 
Since the knowledge of the Hamiltonian spectrum and 
the partition function will be of fundamental importance 
throughout our work and since our goal is not a study 
of full phase diagrams for macroscopic systems, we will 
deal with a simple system composed of two qubits only. 
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Moreover we will work with a simplified XX interaction 
term |22j in order to obtain a particularly easy form for 
Hamiltonian spectrum. 
The Hamiltonian then reads 



(1) 



where A is a coupling constant describing exchange inter- 
actions between the two spins, h is an external magnetic 
field in the z direction, ui {(J2) is the Pauli matrices vec- 
tor for the first (second) atom and a — $1 + $2- 
It is easy to show that this Hamiltonian can be cast in 
the form 



H = -A(S'2 - Si -I)- hS, 



(2) 



of our system whose density matrix p can be cast as fol- 
lows: 



P 



(7) 



Z and H being respectively given by Q and In the 
coupled basis, such an operator reads 
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(8) 



where S ~ ^. The advantage of working with H as given 
by ([2]) is that its diagonalization is now straightforward, 
since this Hamiltonian is diagonal in the coupled ba- 
sis ||5,A/>| = ||1,1),|1,0),|1,-1),|0,0)| of common 

eigenvectors of S*^ and Sz- 
In this ordered basis H reads 



H 



-h 
-A 
h 
A 



(3) 



A. Thermal and entanglement properties, and 
thermal state of the system 

Starting from the knowledge of Hamiltonian eigenval- 
ues it is possible to obtain a closed form for the partition 
function Z of the system which reads 



Z = 2 cosh(/3/i) + 2 cosh(/3A). 



(4) 



Here^=^. 

In view of our target to link thermodynamical quantities 
to entanglement measures, equation Q is of great impor- 
tance since it encompasses all thermodynamical proper- 
ties of our system. We are now able to give an expression 
of the heat capacity. Since (fcs — 1) 



U = - 



d ln(Z) _ h smh{/3h) + A sinh(/3A) 



8/3 



cosh(/3/i) + cosh(/3A) 



(5) 



we obtain 



cv 

,32 



dU 
dl3 



h'-' + X^ + ^iX-hf cosh (^i(A+/i)) + i(A+/i)^ cosh (^(A-h)) 



(cosh(/3/i)+cosh(^(A))2 



(6) 



What we are interested in is the study of a thermal state 



On the other hand we wish to quantify entangle- 
ment between two qubits in our system. In order to 
measure the degree of quantum correlation characteriz- 
ing the two-spin system in the state given by equation 
([7]), it is possible to use many different parameters: one 
of them is the negativity of the state of the system; 
another typical choice is the concurrence. We will use 
this latter parameter thanks to its particularly simple 
expression in terms of Hamiltonian parameters. 
Our entanglement measure will thus be defined as 
C = Max{0, ^'}, where v = ^JJ^_ — ^Jta — ^fj^, — ^fV^- Here 
the /iiS are eigenvalues of the matrix p{ay®<jy)p* (oy^cfy) 
and /^i > /i2 > Ms ^ M4- Evaluating such a quantity, 
after some trivial calculations, we obtain 



C = Max 



{0,;.= |(|sinh^A|-l)}. 



(9) 



We notice that, strictly speaking, the system in its ther- 
mal state is entangled if and only if |A| > In (1 + \/2) = 

A. One could be surprised noticing that such a critical 
value for A does not depend upon h. Nevertheless the 
partition function explicitly depends upon h in such a 
way that, when /i ^ A, C is almost zero even if A > A. 
Only when A > /i the concurrence is significantly dif- 
ferent from zero. Physically this statement means that 
entanglement can arise between the two spins only when 
their mutual interaction is stronger than the magnetic 
one, since a strong magnetic field (with respect to the 
exchange interaction) forces both spins to align with it, 
erasing all correlations between them. 
From equation ^ the concurrence is given in terms of 
three parameters, namely /3, A and h. One way to pro- 
ceed is to look at /i as a parameter characterizing a class 
of functions C/i(/3, A). In this way it is possible to find 
analytical expressions for the functions Aft (C,/3), valid at 
least in a range of values of C. 

What we will indeed be able to do is to find analytical 
expressions for the functions A/i(i^, /3) each of which, as 
long as > 0, coincides with A;i(C,/3). 
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II. HEAT CAPACITY VERSUS 
CONCURRENCE: ANALYTICAL RESULTS AND 
DISCUSSION 



Starting from equation ([9| it is possible to obtain a 
class of inverse functions parameterized by h giving A in 
terms of /3 and v. To this end we write such a set of 
functions as 



fimctions Cy\l3,h') as 



Ch = Max 



{o,< 



cosh ph + + sinh^ 



(|sinh/3A|-l)}. 



(10) 



As we can see from ([9]), and as expected from the form 
of the spectrum (equation ([3|), entanglement is an even 
function of A. This is due to the fact that the two eigen- 
vectors 1 10) and 1 00) associated to the eigenvalues ±A 
manifest the same degree of quantum correlations. 
If we then set x = sinh/3A and a = cosh/3/i and we limit 
our analysis to positive values of A we obtain for v. 



x-1 



a+ 



It is then straightforward to obtain 



(11) 




-^{Xhiiy, P) - hf cosh {I3{\h{y, P) + h)) + 

-{\h{v, P) + hf cosh (/3(A„(z., /3) - h))^ . 

(14) 



Equation ( 14 1 is the main result of this work. It gives 



direct analytical expressions of a thermal parameter in 
terms of an entanglement measure and shows how, once 
the entanglement in the system is known, heat capacity 
has a well defined value. Unfortunately, the converse is 
not true. 

In order to perform an analysis of the link between Cy 
and C we can choose one single function out of the set 
Cy\p^v). In other words, we have to fix the value of 
h with respect to the other parameters involved in ([2|. 
One possible choice is to measure all energies in units 
of h, thus obtaining the function Cy\l3,v). In Fig.l 
we show the behavior of such a function for 9 different 
values of /?. 



av + 
1^^ 



HeatCapacity 



Xhiv.P) = -ArcSinh 

(12) 

As said, as long as we study such a function in the range 
V € [0,1] we actually obtain the exact dependence of A 
upon C. Nevertheless, as shown before, there is a whole 
range of values of A which results in a negative value of v 
or, in other words, which results in a zero concurrence. In 
the whole range A e [0, i In (1 + y/2)] the concurrence as 
a function of A is not an invertible function. To overcome 
such an obstacle we concentrate ourselves on u which, 
when non-negative, coincides with the real concurrence. 
In this sense it is possible to study how heat capacity 
depends upon v in the whole range of values of A, hav- 
ing in mind that the only physically relevant results can 
be obtained by looking at the restriction of the functions 
C'^\l3,v) to the set v e [0,1], while for negative val- 
ues of V nothing can be said about the relation between 
heat capacity and entanglement, since the latter is always 
equal to zero even if the former has not a constant value. 
For fixed temperature and magnetic field the domain of 
the function Cy\l3, v) is then given by 




Inverse Temperatures in units of li 



FIG. 1: Heat Capacity Cv of the system versus Concurrence 
C (all energies are measured in units of magnetic field h) for 
9 difi'erent values of /3: 13 = 0.5 (blue line), /3 = 0.8 (brown 
line), /3 = 1 (red line), /3 = 1.5 (gray line), /3 = 2 (purple 
line), /3 = 3 (orange line), /5 = 4 (black line), /3 = 6 (green 
line) and /3 = 12 (blue dashed line) 



D 



(13) 



as can be easily seen from ^ by evaluating the values of 
1/ for A = and for A — ?► oo. 

By substituting ( 12 ) in ([6| we finally obtain the set of 



From this figure it is possible to notice how heat 
capacity always goes to zero for maximally entangled 
states, irrespectively of temperature. This feature may 
be due to the fact that a strong entanglement prevents 
the system from going towards excited states because of 
the intense correlations between its two parts. 
Another interesting detail worth noticing is that lowering 
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temperature results in a narrower range of negative 
values of v as can be seen from equation (13 1. This 



means that in the limit of very low temperature the 
behavior of heat capacity can be fully described by 
changes in entanglement between the two spins in the 
system. This nice feature is of great interest in view of 
a detailed study of quantum phase transitions in spin 
systems. 

Third, as it is easily seen from Fig.l, heat capacity is 
never a monotonic function of concurrence. Its behavior 
has to be analyzed in detail in order to get an insight 
of the physics behind it. Let us start by looking at the 
low temperature behavior of heat capacity. When the 
state of the system is separable (zero concurrence) heat 
capacity is zero. In Fig. 2 the energy spectrum of our 
system is shown versus A. A zero value of concurrence 
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Hamiltoniaii spectrum 

|1-1> |11> |00> |10> 



Ai = i?ei — Eg becomes much smaller than ksT and 
thus the first excited level is almost as populated as 
the ground state one, in order to raise its energy the 
system has to populate the second excited level. The gap 
between the latter and the ground state is much greater 
than Ai and, at low temperatures, it is also greater 
than fc^r. Thus heat capacity decreases and eventually 
becomes zero. This explains the first maximum in Fig.l. 
When the critical point is crossed, the same reasoning 
applies and thus another maximum appears. Finally 
for very high values of A (corresponding to maximum 
entanglement) the gap becomes wider and wider and 
heat capacity rapidly goes to zero. 

The same kind of argument can be applied to explain 
the high-temperature behavior of heat capacity, starting 
from a small but not zero value (high gap, but not much 
greater than mean thermal energy), growing towards a 
maximum (approaching the critical point) and finally 
decreasing to zero for high values of A. 
Finally, the mid-temperature heat capacity shows a pe- 
culiar behavior for small concurrence (small A): instead 
of increasing with the decreasing gap between ground 
state and first excited levels, heat capacity decreases. 
This can be explained with the same argument employed 
to describe the minimum close to quantum critical 
point: when temperature is not low enough, for small 
A both the ground state and the first excited level are 
significantly populated. This means the system has to 
move to the second excited level in order to raise its 
energy, and the gap between this latter level and the 
ground state increases with increasing A. This explains 
why heat capacity decreases for low concurrence. 



FIG. 2: Hamiltoniaii spectrum versus A. All energies are 
measured in units of h 

can only be obtained for almost zero values of A, when 
the ground state is a factorized one. We see from Fig. 2 
that for very small absolute values of A the gap between 
ground state and all the excited levels is significantly 
greater than the mean thermal energy when for example 
/3 = 12 in Fig.l. Thus the system can not move to an 
excited state, resulting in a zero heat capacity. From 
Fig.l we see how, moving towards greater values of 
concurrence, heat capacity shows two maxima and a 
minimum between them. Such a behavior is typical of 
gapped two level systems and can be seen as related to 
a double Schottky-like anomaly [23l [24] right before and 
right after the quantum critical point. This shows how 
at low temperatures our system actually behaves as a 
two-level one. The two maxima can then be explained 
exploiting Schottky argument, or simply by looking 
at hamiltonian spectrum: raising concurrence means 
raising A (|^ > 0) and thus approaching the quantum 
critical point or, in other words, reducing the gap 
between ground state and first excited level. The system 
can now absorb thermal energy from the environment 
and heat capacity starts growing. Nevertheless, once 



III. CONCURRENCE FROM HEAT CAPACITY: 
AN EXPERIMENTAL PROTOCOL 

It is worth stressing once again that negative values 
of v are associated with a zero concurrence and then for 
V <Q the behavior of heat capacity can not be described 
in terms of entanglement between the two spins, since 
this latter quantity is always zero. 

This means that a measure of heat capacity is not able 
to give a definitive reply to the question about entangle- 
ment between the two particles. Let us focus for instance 
on the curve with /3 = 2 (see Fig. 3). It is possible to iden- 
tify three different ranges in C: the first one is associated 
with a negative value of ^^q^'^^ , meaning that only one 
value of concurrence is associated to each value of heat 
capacity (namely, Cv{v) is invertible in that range). The 
second one is characterized by a bouncing behavior of 
heat capacity and thus Cv{v) is not an invertible func- 
tion. Finally in the third range Cv{v) is again invertible 
and each value of heat capacity is associated to a sin- 
gle value of concurrence. The existence of these ranges 
is a direct consequence of the behavior of the function 
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> 0, 



Cv{X). Indeed, since ^ = ^f^, and since 
the oscillating behavior of Cy {v) is due to the oscillating 
behavior of Cy(A). 

Let us now suppose we actually measure heat capacity for 
a real system described by Hamiltonian M. We suppose 
we can control values of temperature (by for instance 
keeping our system in contact with a reservoir) and of 
external magnetic field. Since our system is fixed, A has 
a constant value which, nevertheless, is not known to the 
experimentalist. We wonder to what extent it might be 
possible to obtain information about microscopic details 
of our system (concurrence, entanglement, values of A) 
by measuring Cy. In order to fix ideas, let us suppose 
we measure all energies in units of magnetic field and 
we fix temperature at /3 = 2. Finally, we measure heat 
capacity to get the value Cmeas2 {Gmeas2 = 0.4 in the ex- 
ample shown in Fig. 3). It is easily seen from Fig.l that 
there exists a situation in which such a value of heat 
capacity can be obtained with three different values of 
concurrence, and by this single measure there is no way 
we can state which one of these three values is the real 
one. In other words there exist three different systems, 
physically characterized by different values of A (since h 
and P are fixed), showing the same heat capacity. The 
possible values of A compatible with the measured value 
Gmea82 ^rc casily obtainable from the knowledge of the 
function Cy(A, /3 = 2) (l6]) or from the knowledge of both 
Cy(i/,/3 = 2) PI and A(z^, /3 = 2) (pi). 
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FIG. 3: Heat Capacity Cv of the system versus parameter 
V (all energies are measured in units of magnetic field h) for 
2 different values of /3: /3 = 2 (blue line) and /3 = 3 (red 
line). In figure are shown by dashed lines the changes in heat 
capacity and concurrence when going from /3 = 2 to /3 = 3 
for the three possible values of A associated to Cmeas2 ~ 0.4. 
It can be noticed from this figure how the changes in heat 
capacity strongly depend on the starting point on the curve 
with (5 — 2. Heat capacity gaps Cv[P = 2) — Cv{P ~ 3) are 
sensibly difi'erent for the three possible values of A discussed in 
the text, enabling us to distinguish between the three possible 
physical systems associated to them 

We thus have three possible systems respectively char- 
acterized by values Ai, A2 and A3. Looking at the ex- 
ample of Fig. 3 with Cmeas2 = 0.4 these three values are 



Ai = 0.76054, A2 = 1.59078 and A3 = 2.68228. The ques- 
tion arises whether it may be possible to identify which 
one of these three values is the real one. 
Let us now suppose we measure heat capacity for the 
same system having a different temperature, say /3 = 3. 
Since changing temperature would not affect the real 
value of A characterizing our physical system, and since 
this value belongs to the set Sa = {Ai, A2, A3}, we would 
for sure measure one of the three values Cy (Ai, (3 = 3) — 
Ci, Cy(A2,/3 = 3) = C2 or Cy(A3,/3 = 3) = C3. Since 
these values are in general very different from each other 
it is possible for us to identify which one is the one we 
actually measured (Cmeass) in our second measuring pro- 
cess. Indeed there would be a set St — {A4,A5,A6} of 
values of A associated with the measured values Cmeas3 
of Cy. Following Fig. 3 we might suppose the new value 
of heat capacity to be Cmeass = 0.164 and the set 
St = {0.853576,1.23367,2.682}. There wiU be at least 
one common value between Sa and Sb- It is then possible 
in this way to identify the real value Xp of the coupling 
constant A, since it will be the one for which both of the 
following equalities would hold 



Cy (A — Ap, /3 — 2) — Cmeas2 
Cy (A = Ap, (5 — 3) — Cmeas3- 



(15) 



In the example exposed above the only common value 
between the sets Sa and Sb is found to be Ap = 2.682. 
It is of course possible that more than one of the values 
{Ai, A2, A3} is associated with the measured heat capac- 
ity Cmeas3- In this casc it would be enough to measure 
Cy for a third temperature, to find the new set (Sc) of 
possible values of A associated with this new value of 
heat capacity and to look for the common element of Sa, 
Sb and Sc- This procedure can be employed until the 
physical value of A is uniquely identified. 
Once we know the actual value Ap of A we know every- 
thing about our system. In particular it is possible to 
calculate concurrence (and thus the entanglement in our 
system) using ^ which, as shown, is invertible and thus 
gives a unique value for C. For the situation shown in 
Fig. 3, exploiting the value Ap — 2.682, the concurrence 
is found to be C = 0.956905. 

With simple measurements of heat capacity of the 
two-qubit system we are thus able to measure the 
entanglement as given by concurrence and any other 
microscopic property of the system. 



IV. CONCLUSIVE REMARKS 

It might be interesting to compare our results about 
heat capacity and entanglement with the one shown in 
|15) , where an upper bound for heat capacity of systems 
in entangled states has been obtained. It is shown that 
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there exists a maximal value of heat capacity for entan- 
gled states, such that any system showing heat capacity 
greater than this value is for sure in a disentangled state. 
Briefly, the so-called separable bound on heat capacity is 
a direct consequence of the existence of a separable bound 
on internal energy. Indeed it has been shown that for 
composite systems every separable state is associated to 
an energy greater than a certain value Eb- Moreover an 
explicit form for this separable bound on internal energy 
has been given for magnetic spin systems [25]. Applied 
to our system such a bound becomes, for any separable 
state, U > Eb with 



1 -\h\ 



^ \h\<2\X\ 



\h\>2\X\. 



(16) 



Exploiting the result (obtained in [TS]) giving the sepa- 
rable bound on heat capacity for gapped systems (valid 
for low temperature only) as 



Eg) - Cy 



(17) 



where Eg is the ground state energy, A is the energy gap 
between ground state and the first excited level and Eb 



is given by (16), we can obtain an expression for such 



a bound valid for the low temperature behavior of the 
system analyzed in this work, which reads 



h), 



= < 



(-A-M(-|A|-^ 
(A + ;.)(-|A|-M^ + 
{\ + h){-\h\ + h), 
i-\ + h){-\h\ + h), 

(-A + M(-|A|-^ + /i), 
(A-/^)(-|A|-g^ + A), 



A), \<~h 



-h < X < 



< A < 



< A < f 

1 < A < /i 
X>h 



(18) 



Notice that Cy as a function of A and h is always nega- 
tive, while by definition heat capacity is a positive quan- 
tity. This is a consequence of the fact that energy sep- 
arable bound is for our system always smaller than the 
ground state energy. Thus this method would give no 
definite reply to the question about entanglement at low 
temperature since (17) always holds but its violation is a 



sufficient but not necessary condition for the existence 
of quantum correlations. This shows how this latter 
method and the one reported in this paper, both based 
on measures of heat capacity, are not equivalent. It is 
worth noticing that, while the method based on sepa- 
rable bounds is (when possible) only able to detect the 
presence of entanglement, the one shown in this work is 
able to measure entanglement by means of an (indirect) 
measure of concurrence. 

It is worth stressing that the low temperature behav- 
ior of entanglement is the one characterizing the criti- 
cal behavior of the system and is thus of great interest; 
moreover when the temperature is high enough, thermal 
entanglement is usually destroyed. Thus the comparison 
presented above, even if valid only for low temperature, 
is yet full of physical meaning. 

Acknowledgements 



One of the authors (H.N.) acknowledges the support 
from JSPS, Japan through a Grant-in- Aid for Scientific 
Research (C). 



References 



[1] 


A. Osterloh, L. Amico, G. Falci, R. Fazio, Nature 410, 




250404 (2004) 


608 (2002) 


[11] 


G. Mahler, J. Gemmer, M. Michel, Physica E 29, 53 


[2] 


T.J. Osborne, M.A. Nielsen, Phys. Rev. A 66, 032110 




(2005) 




(2002) 


[12] 


M. Hartmann, G. Mahler, O. Hess, Phys. Rev. Lett. 93, 


[3] 


M. Greiner, 0. Mandel, T. Esslinger, T. W. Hansch, I. 




080402 (2004) 


Bloch, Nature 415, 39 (2002) 


[13] 


J.-Q. Liao, J.-F. Huang, L.-M. Kuang, arXiv: 1007.4296 /I 


[4] 


S. Sachdev, Quantum Phase Transitions, Cambridge Uni- 




(2010) 




versity Press, 2008 


[14] 


M. C. Bafiuls, J. L Cirac, M. B. Hastings, Phys. Rev. 


[5] 


H.-M. Kwok, C.-S. Ho, S.-J. Gu, Pliys. Rev. A 78, 062302 




Lett 106, 050405 (2011) 




(2008) 


[15] 


M. Wiesniak, V. Vedral, C. Brukner, Phys. Rev. B 78, 


[6] 


W. Son, L. Amico, F. Plastina, V. Vedral, Phys. Rev. A 


064108 (2008) 




79, 022302 (2009) 


[16] 


R. Brustein, A. Yarom, Phys. Rev. D 69, 064013 (2004) 


[7] 


J. Vidal, G. Palacios, R. Mosseri, Phys. Rev. A 69, 


[17] 


M. Asoudeh, V. Karimipour, Phys. Rev. A 70, 052307 




022107 (2004) 


(2004) 


[8] 


M. C. Arnesen, S. Bose, V. Vedral, Phys. Rev. Lett. 87, 


[18] 


H. T. Quan, Z. D. Wang, C. P. Sun, Phys. Rev. A 76, 




017901 (2001) 


012104 (2007) 


[9] 


S.-J. Gu, G.-S. Tian, H.-Q. Lin, Phys. Rev. A 71, 052322 


[19] 


J. Ren, S. Zhu, X. Hao, J. of Phys. B: At. Mol. Opt. 




(2005) 


Phys. 42, 015504 (2009) 


[10] 


L.A. Wu, M.S. Sarandy, D.A. Lidar, Phys. Rev. Lett. 93, 


[20] 


S. CampbeU, M. Paternostro, Phys. Rev. A 82, 042324 



7 



(2010) 

[21] P. Hanggi, G.-L. Ingold, P. Talkncr, New J. of Phys. 10, 

115008 (2008) 
[22] X. Wang, Phys. Rev. A 64, 012313 (2001) 
[23] W. Greiner, L. Neise, H. Stocker, Thermodynamics and 

statistical mechanics, Springer, 1997, p. 221 
[24] A Schottky anomaly (also known as Schottky heat capac- 



ity) is the behavior of heat capacity versus temperature, 
typical of a two level system, characterized by a central 
maximum followed by an asymptotic decreasing toward 
zero; the only difference here is that temperature is fixed 
and the only parameter being varied is A. 
[25] G. Toth, Phys. Rev. A 71, 010301(R) (2005) 



